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PREFACE 


The  study  described  in  this  report  was  conducted  as  part  of  the 
Environmental  and  Water  Quality  Operational  Studies  (EWQOS)  Program, 

Work  Units  31594  (IB.l)  and  31595  (IC.l).  The  EWQOS  Program  is  spon¬ 
sored  by  the  Office,  Chief  of  Engineers  (OCE) ,  US  Army,  and  is  assigned 
to  the  US  Army  Engineer  Waterways  Experiment  Station  (WES)  under  the 
purview  of  the  Environmental  Laboratory  (EL).  The  uCE  Technical  Moni¬ 
tors  for  EWQOS  were  Mr.  Earl  Eiker,  Dr.  John  Bushman,  and  Mr.  James  L. 
Gottesman. 

This  is  an  interim  report  dealing  with  work  conducted  to  December 
1982.  Work  on  evaluation  of  the  model  will  continue  through  Fiscal  Year 
1984. 

The  work  reported  herein  was  conducted  by  Drs.  Joseph  H.  Wlosinski 
and  Carol  D.  Collins  of  the  Water  Quality  Modeling  Group  (WQMG) ,  Ecosys¬ 
tem  Research  and  Simulation  Division  (ERSD) ,  EL.  The  draft  report  was 
reviewed  by  Dr.  James  L.  Martin,  Dr.  Stephen  P.  Schreiner,  and 
Mr.  Mark  S.  Dortch,  WQMG.  The  study  was  conducted  under  the  direct 
supervision  of  Mr.  Don  L.  Robey,  Chief,  ERSD,  and  under  the  general 
supervision  of  Dr.  John  Harrison,  Chief,  F,L.  Program  Manager  of  EWQOS 
was  Dr.  J.  L.  Mahloch,  EL.  Editorial  review  was  performed  by 
Ms.  Jessica  S.  Ruff  of  the  WES  Publications  and  Graphic  Arts  Division. 

Director  of  WES  was  COL  Allen  F.  Grum,  USA.  Technical  Director 
was  Dr.  Robert  W.  Whalin. 

This  report  should  be  cited  as  follows: 

Wlosinski,  J.  H. ,  and  Collins,  C.  D.  1985.  "Analysis  and  Revi¬ 
sion  of  a  Reservoir  Water  Quality  Model,"  Technical  Report 

E-85-13,  US  Army  Engineer  Waterways  Experiment  Station,  Vicksburg, 

Miss . 
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ANALYSIS  AND  REVISION  OF  A  RESERVOIR 


WATER  QUALITY  MODEL 


PART  I:  INTRODUCTION 
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1.  CE-QUAL-R1  is  a  one-dimensional  model  that  is  being  developed 
by  the  Corps  of  Engineers  to  realistically  predict  and  assess  the  ef¬ 
fects  of  engineering  activities  on  reservoir  water  quality.  This  report 
deals  with  the  analysis  and  revision  of  that  model  up  to  December  1982. 
Further  evaluation  is  scheduled  to  continue  through  Fiscal  Year  1984. 

The  purposes  of  this  work  were  to  ensure  that  the  coding  of  the  model  is 
correct  and  that  model  predictions  are  suitable  for  the  needs  of  Corps 
of  Engineers  District  and  Division  offices. 

2.  The  evaluation  of  model  predictions  from  CE-QUAL-R1  was  sepa¬ 
rated  into  three  tasks,  with  each  task  being  investigated  by  a  separate 
group.  The  three  tasks  dealt  with  the  predictions  of  temperature  (John¬ 
son  and  Ford  1981),  anaerobic  materials  (Zimmerman  1985),  and  other 
chemical  and  biological  variables.  In  addition,  an  evaluation  of  the 
model  for  use  by  the  Aquatic  Plant  Control  Research  Program  was  reported 
by  Wlosinski  (1981).  Evaluation  of  model  predictions  is  also  being  car¬ 
ried  out  using  data  collected  on  a  reservoir  with  different  morphometric 
and  biological  attributes  (Wlosinski  and  Collins  1985).  Variables  exam¬ 
ined  in  the  analysis  and  reported  here  were  dissolved  organic  matter 
(DOM),  algae,  zooplankton,  orthophosphate-phosphorus  (P04-P) ,  ammonia- 
nitrogen  (NH4-N) ,  nitrite  plus  nitrate-nitrogen  (N02-N  +  N03-N) ,  total 
inorganic  carbon,  oxygen,  pH,  alkalinity,  and  total  dissolved  solids 
(TDS) .  Data  for  the  analysis  were  collected  at  DeGray  Lake,  Arkansas, 
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in  1979  and  1980. 


PART  II:  BACKGROUND  INFORMATION 


Model  Description 

3.  CE-QUAL-R1  is  a  one-dimensional  numerical  reservoir  model  that 
simulates  water  quality  variables  in  the  vertical  direction.  The  thick¬ 
ness  of  each  horizontal  layer  is  dependent  upon  the  balance  of  inflowing 
and  outflowing  water,  which  permits  accurate  mass  balancing  and  reduces 
numerical  dispersion  during  periods  of  large  inflow  and  outflow. 

4.  Inflowing  waters  are  distributed  vertically  based  on  density 
differences,  which  allows  simulation  of  surface  flows,  interflows,  and 
underflows.  Water  density  is  dependent  on  temperature  and  the  concen¬ 
trations  of  dissolved  and  suspended  solids.  Outflowing  water  is  with¬ 
drawn  from  layers,  considering  density  stratification,  using  the  selec¬ 
tive  withdrawal  algorithms  of  Bohan  and  Grace  (1973).  Reservoir 
outflows  can  be  specified  by  using  operation  records  or  the  user  may  opt 
to  have  the  model  choose  flows  from  ports  in  order  to  match  a  downstream 
target  temperature. 

5.  The  heat  budget  includes  the  components  of  short-  and  long¬ 
wave  radiation,  back  radiation,  reflected  solar  and  atmospheric  radia¬ 
tion,  evaporative  loss,  conductive  heat  transfers,  and  gain  or  loss 
through  inflows  and  outflows.  Vertical  transport  of  thermal  energy  and 
mass  is  achieved  through  entrainment  and  turbulent  diffusion.  Entrain¬ 
ment  determines  the  depth  of  the  upper  mixed  layer  and  the  onset  of 
stratification.  It  is  calculated  from  the  turbulent  kinetic  energy 
influx  generated  by  wind  shear  and  convective  mixing  using  an  integral 
energy  approach  (Johnson  and  Ford  1981).  Turbulent  diffusion  is  a  two- 
way  transport  process  that  is  incorporated  using  a  turbulent  or  eddy 
diffusion  coefficient,  which  is  dependent  on  the  wind  speed,  magnitude 
of  inflows  and  outflows,  and  density  stratification. 

6.  The  prediction  of  water  quality  is  based  upon  simulation  of 
the  interaction  of  numerous  biological  and  chemical  constituents. 

Forces  that  directly  affect  the  simulation  of  the  biological  and  chemi¬ 
cal  constituents  are  temperature,  irradiation,  wind  speed,  inflow  and 


outflow  rates,  and  Inflowing  and  outflowing  masses.  The  physical  dis¬ 
tribution  of  mass  is  dependent  upon  the  diffusive  and  convective  pro¬ 
cesses  described  above  and  on  settling  processes.  Besides  the  physical 
processes,  water  quality  variables  in  the  model  can  be  affected  by  pho¬ 
tosynthesis,  respiration,  decay  or  decomposition,  ingestion,  egestion, 
nonpredatory  mortality,  and  harvest.  Table  1  lists  the  process  interac¬ 
tions  between  variables  which  were  included  in  CE-QUAL-R1  at  the  begin¬ 
ning  of  this  study.  A  definitive  description  of  the  model  can  be  found 
in  a  User's  Manual  (Environmental  Laboratory  1982b). 

DeGray  Lake 

7.  DeGray  Lake  is  located  in  the  Ouachita  Mountains  in  south- 

central  Arkansas.  Reservoir  length  is  32  km,  with  a  maximum  depth  of 

8  3 

57  m.  Volume  at  normal  pool  is  7.91  x  10  m  ,  with  a  surface  area  of 
7  2 

5.34  x  10  m  .  DeGray  Lake  is  dendritic,  with  a  shoreline  development 
index  of  12.8.  Project  purposes  include  flood  protection,  hydropower 
with  pumped  storage,  recreation,  water  supply,  and  low-flow  augmentation 
of  the  Ouachita  River. 

8.  The  reservoir  is  formed  by  the  DeGray  Dam,  which  is  located 
12.7  km  upstream  of  the  mouth  of  the  Caddo  River.  The  dam  is  an  earth- 
fill  structure  with  a  crest  elevation  of  138  m  above  mean  sea  level. 
Water  is  withdrawn  through  a  multilevel  outlet  structure.  From  first 
filling  in  1971  to  March  1979,  water  was  withdrawn  from  the  surface, 

after  which  time  the  withdrawal  level  was  lowered  12  m.  The  watershed 

2 

above  the  dam  is  1,170  km  ,  of  which  approximately  69  percent  is  clas¬ 
sified  as  forested,  30  percent  agricultural,  and  only  1  percent  urban 
(Perrier,  Harris,  and  Ford  1977). 

Evaluation  Procedures 


9.  Procedures  for  evaluating  CE-QUAL-Rl  are  included  in  a  report 
by  Wlosinski  (1984).  Methods  generally  follow  recommendations  made  dur¬ 
ing  a  workshop  on  verification  of  water  quality  models  which  was 


convened  by  the  US  Environmental  Protection  Agency  (1980).  Participants 
of  the  workshop  recommended  evaluation  of  both  software  and  model  pre¬ 
dictions.  Although  they  encouraged  the  use  of  statistical  techniques, 
they  did  not  recommend  any  statistics  nor  did  they  believe  that  statis¬ 
tical  techniques  should  supersede  engineering  judgment. 

Evaluation  of  software 

10.  The  following  methods  were  used  to  evaluate  software: 

a.  Model  equations  were  checked  for  dimensionality. 

b.  Model  predictions  were  checked  for  stability. 

c.  Model  output  was  checked  for  conservation  of  mass. 

d.  All  variables  were  checked  to  ensure  that  they  were  ini¬ 
tially  set. 

e.  Predictions  were  checked  using  different  time  steps. 

ff.  As  suggested  by  Mihram  (1972)  and  Lawler  (1980),  model 
predictions  were  checked  to  make  sure  they  were  reason¬ 
able  after  changing  values  of  driving  variables  (boundary 
conditions) . 

g.  Problems  reported  by  others  using  the  model  were  thor¬ 
oughly  investigated. 

11.  The  majority  of  these  tests  were  performed  in  1980,  with  the 
corrections  being  made  to  the  code  prior  to  the  release  of  the  original 
User's  Manual  (Environmental  Laboratory  1982a). 

Evaluation  of  model  predictions 

12.  Model  predictions  were  evaluated  both  graphically  and  statis¬ 
tically  using  data  from  DeGray  Lake  for  1979  and  1980.  Although  a  num¬ 
ber  of  statistics1  tests  were  available  (see  K'losinski  1984),  most  of 
the  comparisons  for  the  study  were  made  using  the  Reliability  Index  (RI) 
of  Leggett  and  William-  t  il,8I)  .  As  the  predicted  versus  measured  values 
diverge,  the  Rl  becomes  larger.  because  the  RI  does  not  depend  on 
whether  the  observed  or  predicted  value  is  greater,  and  since  it  is 
scale  variant,  it  appears  to  be  the  best  statistic  for  aggregating  and 
comparing  results  of  different  variables. 

13.  Statistics  were  calculated  for  vertical  profiles  for  each  var¬ 
iable  on  those  dates  when  data  were  collected,  and  for  each  variable 
summed  over  depths  and  dates.  The  value  of  the  Rl  is  1.0  in  the  case  of 
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perfect  prediction.  If  all  comparisons  were  within  a  factor  of  two  of 
each  other,  the  R1  would  be  2.0.  An  RI  of  10.0  signifies  that  the  dif¬ 
ferences  between  measured  and  predicted  values  were  an  average  of  one 
order  of  magnitude  apart;  an  RI  of  100.0  signifies  the  values  were  two 
orders  of  magnitude  apart.  The  RI  was  used  for  comparing  predictions 
for  different  model  algorithms,  processes,  structures,  or  calibration 
values.  These  comparisons  were  made  using  data  collected  at  Station  4, 
located  near  the  dam.  This  was  done  because  model  predictions,  in  gen¬ 
eral,  better  represent  conditions  in  the  main  pool  of  the  reservoir,  and 
more  data  were  collected  at  that  station. 

14.  Although  the  model  is  one-dimensional  and  considers  each 
layer  as  a  completely  mixed  reactor,  the  reader  should  be  aware  that 
this  is  not  the  case  in  the  real  system.  Nutrient  inputs  to  a  modeled 
reservoir  layer  become  instantaneously  mixed  throughout  that  layer, 
whereas  in  the  real  system  a  continuum  of  changing  conditions  can  be 
seen  as  one  progresses  from  the  headwater  area  to  the  dam  (Thornton 
et  al.  1980) . 

19.  In  addition  to  comparing  the  mass  of  predicted  and  measured 
state  variables,  fluxes  predicted  by  CE-QUAI.-R1  were  also  compared  with 
measured  values.  This  was  done  in  order  to  ensure  a  more  reliable 
model,  for  it  is  possible  to  predict  the  s"  concentrations  with  dif¬ 
ferent  sets  of  coefficients  (Wlosinski  1985). 


PART  III:  EVALUATION  OF  SOFTWARE 


Stabilitv  of  Model  Predictions 


16.  When  the  model  predictions  were  originally  checked  for  sta¬ 
bility,  problems  occurred  for  those  variables  associated  with  the  bottom 
of  the  reservoir.  At  the  time  this  problem  was  investigated,  these  var¬ 
iables  were  organic  sediment  and  benthos;  however,  since  then,  a  number 
of  variables  representing  anaerobic  materials  have  been  added.  Sample 
results  obtained  from  the  original  model  for  the  bottom  three  layers  of 
sediment  are  shown  in  the  upper  part  of  Table  2.  Although  the  value  for 

sediment  should  be  fairly  stable  over  short  periods  of  time,  predicted 
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values  for  the  bottom  three  layers  ranged  from  0  to  562  g/irT  for  the 
first  4  days  of  simulation.  The  solution  scheme  was  found  to  be  the 
cause  of  this  problem. 

17.  Although  initial  conditions  for  benthos  and  sediment  were  re¬ 
ported  as  grams  per  square  meter,  those  values  were  changed  to  a  concen¬ 
tration  (grams  per  cubic  meter)  by  "dispersing"  the  variables  into  the 
water  of  t tie  adjacent  layer.  Since  the  ratio  of  the  mixture  changed, 
the  concentration  of  the  variable  changed,  becoming  much  greater  in  the 
deeper  layers.  An  implicit  integration  scheme  was  used  to  solve  for  new 
concentrations.  To  solve  the  stability  problem,  the  units  for  the  vari¬ 
ables  were  changed  from  grams  per  cubic  meter  to  grams  per  layer,  and 
the  1'iiler  technique  was  used  to  solve  for  new  values.  Predictions  of 

•  1  i  -  on:  :  or  da  vs  in  the  bottom  1  avers  for  the  imr roved  model  ranged 

•  •twee.  ’  ■  and  li.l  g.'m  and  were  considered  satisfactory.  Results 
after  corrections  to  the  model  are  listed  in  the  bottom  portion  of 


Conservation  of  Mass 

iH.  A  thorough  test  for  conservation  of  mass  was  performed  using 
data  from  lake  fonv.iv,  Florida  (Wlosinskf  l^Ri).  For  a  conservative 
substance  in  the  model,  such  as  total  dissolved  solids,  the  results  were 
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satisfactory  after  a  1-year  simulation.  The  balance,  which  is  a  compar¬ 
ison  of  the  initial  mass  plus  all  positive  fluxes  of  the  reservoir 
boundary  with  the  final  mass  plus  all  negative  fluxes,  was  within 
1.7  percent.  For  nonconservative  substances,  such  as  phosphorus,  the 
initial  results  were  not  satisfactory.  In  a  series  of  three  calibration 
simulations,  the  balance  for  phosphorus  was  between  4.5  and  147.5  per¬ 
cent  (Table  3).  The  three  calibration  simulations  used  exactly  the  same 
data  set  except  for  the  estimates  for  three  coefficients.  The  algae 
half-saturation  coefficient  for  phosphorus  was  varied  from  0.005  to 
0.003;  the  algal  settling  rate,  from  0.4  to  0.2;  and  zooplankton  assimi¬ 
lation,  from  0.33  to  0.27.  The  range  of  these  estimates  is  reasonable 
and  can  be  found  in  Jorgensen  (1979). 

19.  The  imbalance  of  mass  was  caused  by  three  factors.  First, 
the  ratios  of  carbon,  nitrogen,  and  phosphorus  were  different  for 
different  state  variables.  Although  this  difference  occurs  in  nature, 
the  model  is  not  complex  enough  to  represent  all  of  the  processes  in 
nature  needed  to  account  for  an  exact  mass  balance.  For  example,  the 
ratios  of  carbon,  nitrogen,  and  phosphorus  for  zooplankton  and  algae  are 
usually  different.  When  zooplankton  eat  algae,  they  do  not  take  on  the 
carbon,  nitrogen,  and  phosphorus  ratios  of  algae,  but  instead  excrete 
these  elements  in  different  ratios.  In  the  original  model,  when 
zooplankton  ate  algae,  elements  were  created  or  destroyed  since  algae 
and  zooplankton  had  different  carbon,  nitrogen,  and  phosphorus  ratios. 

20.  Three  alternatives  were  available  to  correct  this  problem: 

(a)  the  model  could  have  been  made  more  complex  to  better  represent 
natural  processes  dealing  with  the  ratio  of  elements,  (b)  the  assumption 
could  be  made  that  different  organic  state  variables  have  the  same 
stoichiometric  ratio,  or  (c)  the  assumption  could  be  made  that  the  prob¬ 
lem  did  not  cause  significant  errors.  Because  the  first  alternative 
would  have  resulted  in  a  model  too  complex  for  the  needs  of  the  Corps  of 
Engineers,  and  because  the  mass  balance  was  considered  significant,  the 
second  alternative  was  chosen.  Thus,  a  simplified  assumption  has  been 
made  in  the  model  that  a  constant  stoichiometric  ratio  exists  between 


the  elements  of  carbon,  nitrogen,  and  phosphorus  for  different  organic 
variables . 

21.  The  second  problem  causing  the  mass  imbalance  was  due  to  the 
solution  scheme.  The  model  solves  the  equations  dealing  with  each  state 
variable  in  a  sequential,  manner.  In  effect,  this  scheme  solves  coupled 
differential  equations  in  an  uncoupled  manner.  The  flux  between  com¬ 
partments  is  actually  solved  twice:  once  for  computations  for  the  donor 
compartment  and  once  for  the  receiving  compartment.  The  flux  computed 
by  the  model  can  actually  be  different  for  the  two  compartments  since  a 
new,  updated  value  for  the  mass  of  the  donor  compartment  was  used  in 
calculating  the  flux  in  the  receiving  compartment.  This  problem  was  not 
found  to  be  a  major  factor  causing  mass  imbalance. 

22.  The  third  problem  caused  the  majority  of  the  mass  imbalance. 
Under  certain  conditions,  more  material  was  predicted  to  leave  a  com¬ 
partment  than  was  contained  in  that  compartment.  On  those  occasions, 
the  model  arbitrarily  changed  the  predicted  negative  concentration  to 
either  zero  or  a  small  positive  concentration.  This,  in  effect,  created 
mass  that  was  then  used  by  other  compartments  in  the  model,  just  as  if 
the  addition  entered  the  reservoir  along  with  the  upstream  flow.  During 
the  Lake  Conway  application,  these  arbitrary  additions,  termed  the  nega¬ 
tive  hedge,  were  totaled  in  order  to  assess  their  significance 

(Table  3).  As  can  be  seen  for  simulation  3  (original  model),  the  amount 
of  phosphorus  added  by  way  of  the  negative  hedge  was  more  than  the  total 
initial  mass  or  the  total  amount  in  the  inflow. 

23.  The  algorithms  for  those  processes  which  caused  the  majority 
of  the  negative  hedge  problem  were  altered  so  that  the  maximum  amount  of 
material  that  could  leave  a  compartment  was  the  amount  present  at  the 
beginning  of  that  time  step.  The  balance  for  the  improved  model  was 

3.1  percent,  which  compared  with  147.5  percent  when  using  the  original 
model  with  the  third  data  set.  The  value  of  mass  added  by  the  negative 
hedge  fell  from  120,000  kg  to  only  20  kg.  The  amount  of  phosphorus 
added  by  way  of  the  negative  hedge  and  the  balance  percentage  were  even 
lower  when  using  the  coefficients  from  simulations  1  and  2  with  the 
improved  model. 


Other  Tests 


24.  Other  tests  performed  on  CE-QUAL-R1  included  checking  the 
dimensions  of  each  equation,  comparing  predictions  using  different  time 
steps,  changing  driving  variables  while  making  sure  predictions  were 
reasonable,  and  checking  to  make  sure  that  initial  values  were  supplied 
for  all  model  variables.  Each  of  these  tests  pinpointed  errors  within 
the  model.  The  majority  of  errors  were  easily  corrected.  Corrections 
ranged  from  changing  misspelled  variable  names  to  replacing  algorithms. 
The  more  significant  changes  included: 

a.  Replacing  an  approximate  solution  for  calculating  eleva¬ 
tion  as  a  function  of  volume  with  an  exact  solution. 

b.  Correcting  the  algorithm  for  estimating  per  hour  rates 
from  per  day  rates  used  in  the  Euler  solution  scheme. 

c.  Changing  the  algorithm  for  interpolating  initial 
conditions. 

d.  Correcting  algorithms  that  dealt  with  the  settling  of 
algae  and  detritus. 

e.  Making  density  a  function  of  dissolved  and  suspended 
solids  in  addition  to  temperature. 


PART  IV:  EVALUATION  OF  MODEL  PREDICTIONS 


Original  Model 

25.  Work  reported  in  this  section  of  the  report  commenced  after 
improvements  from  PART  III  were  made  to  the  model  and  after  satisfactory 
thermal  profiles  were  predicted  for  DeGray  Lake  (Johnson  and  Ford  1981). 
The  work  of  Johnson  and  Ford  provided  temperature  predictions  (from  CE- 
THERM-R1,  the  thermal  analysis  section  of  CE-QUAL-R1)  that  had  an  aver¬ 
age  RI  of  1.11. 

26.  Variables  for  which  data  were  available  from  DeGray  in  1979 
and  1980  were  algae,  zooplankton,  dissolved  organic  matter,  POA-P, 

NHA-N,  N02-N  +  N03-N,  oxygen,  pH,  alkalinity,  and  TDS.  In  addition, 
data  for  inorganic  carbon  were  available  for  1979.  Two  different  esti¬ 
mates  for  algae  were  available  for  both  years.  One  estimate  was  based 
on  dry  weight  values  in  which  zooplankton  were  handpicked  from  a  grab 
sample  with  the  sample  then  being  dried.  The  other  estimate  was  based 

on  chlorophyll  a  data.  For  1979  and  1980,  the  dry  weight  data  averaged 
~  3 

2.27  and  1.68  g/m  ,  respectively.  Values  estimated  from  the  chlorophyll 
data,  using  the  average  concentration  and  a  conversion  value  from 

3 

Spangler  (1969),  were  0.56  and  0.29  g/m  .  Although  the  dry  weight  data 
were  used  initially,  estimates  for  the  dry  weight  may  have  been  high  due 
to  detritus  incorporated  in  the  samples.*  For  this  reason,  and  since 
more  chlorophyll  a  data  were  available,  final  comparisons  were  made 
using  the  chlorophyll  a  estimates. 

27.  Results  for  1979  obtained  from  the  original  model,  with  cor¬ 
rections  as  described  in  Part  III,  are  presented  in  Table  A  and  Fig¬ 
ure  1.  These  results  were  obtained  after  a  number  of  calibration  simu¬ 
lations  were  completed.  The  graphs  in  Figure  1  also  contain  the  results 
after  all  improvements  were  made  to  the  model  (final  model).  (These 


improvements  are  described  later  in  this  section.)  In  those  cases  where 
the  predictions  from  the  original  and  final  models  were  similar,  predic¬ 
tions  from  only  one  are  presented.  If  the  predictions  are  not  similar, 
the  solid  line  represents  predictions  from  the  final  model  and  the 
dashed  line  represents  predictions  from  the  original  model.  Although 
graphs  were  obtained  for  all  variables  for  a  number  of  days  when  mea¬ 
sured  data  were  available,  only  those  graphs  viewed  at  best  and  worst 
cases  (as  determined  by  the  RI)  are  presented.  The  observed  and  pre¬ 
dicted  means  as  well  as  the  RI  are  presented  in  Table  A.  The  statistics 
presented  were  calculated  over  all  depths  and  time  periods. 

28.  Examination  of  results,  from  the  original  model,  showed  that 
the  predicted  POA-P  and  NHA-N  values  were  not  at  all  satisfactory,  espe¬ 
cially  in  the  hypolimnion  during  the  latter  months  of  the  year.  This, 
in  turn,  caused  high  predictions  for  algae  and  inorganic  carbon  in  the 
epilimnion  as  the  mixed  layer  became  deeper  in  the  fall.  This  problem 
persisted  through  all  calibration  simulations,  indicating  that  the  model 
structure  probably  was  not  correct.  Results  obtained  for  the  1980  con¬ 
firmation  data  set  also  showed  high  predictions  for  phosphorus  and  ni¬ 
trogen.  Results  for  the  confirmation  data  sets  are  presented  in  Fig¬ 
ure  2  and  Table  5.  The  format  is  the  same  as  for  the  1979  data  set. 

29.  Further  attempts  to  calibrate  the  model  did  not  produce  bet¬ 
ter  results.  At  this  point,  the  decision  was  made  to  change  the  struc¬ 
ture  and  algorithms  of  the  model  to  see  if  results  could  be  improved. 

Any  changes  that  simplified  the  model  would  be  accepted  even  if  the  re¬ 
sults  were  not  significantly  better.  Changes  that  made  the  model  more 
complex,  and  therefore  harder  to  use,  would  only  be  accepted  if  they 
produced  significantly  better  results. 


Collapsing  Compartments 

30.  The  first  changes  that  were  made  attempted  to  simplify  the 
model.  Changes  were  made  one  at  a  time,  after  which  graphic  and  statis¬ 
tical  analyses  were  performed  with  the  results  of  the  1979  data  set. 


31.  Three  changes  were  made  that  simplified  the  model: 

a.  Three  fish  compartments  were  united  to  form  one 
compartment. 

b.  The  benthos  and  sediment  compartments  were  united. 

c.  Nitrite  and  nitrate-nitrogen  compartments  were  combined 
to  form  one  compartment. 

32.  Diagrams  of  the  new  model  structure  which  incorporated  these 
changes  are  presented  as  Figures  3-5.  Graphic  comparisons  of  all  vari¬ 
ables  after  each  change  appear  to  be  nearly  identical  to  the  graphs 
after  final  calibration.  Only  slight  differences,  before  and  after  com¬ 
bining  variables,  were  noted  in  the  RI ,  and  the  overall  average  of  3.43 
for  all  variables  was  maintained.  All  coefficients  remained  the  same 
except  those  directly  involved  in  combining  variables.  After  all 
changes  were  made,  the  1980  confirmation  data  set  was  again  used  to  sim¬ 
ulate  DeGray  Lake.  As  happened  for  the  1979  data  set,  graphs  for  all 
variables  appear  to  be  nearly  identical  when  compared  with  the  confirma¬ 
tion  simulation  after  final  calibration.  The  overall  RI  was  slightly 
worse,  going  from  3.49  to  3.56.  Because  the  predictions  were  nearly  the 
same  before  and  after  the  changes,  and  since  the  simplified  version  used 
35  fewer  coefficients,  the  changes  were  made  permanent. 

Adsorption  of  Phosphorus  and  Nitrogen 

33.  The  process  of  adsorption,  by  which  materials  in  the  soluble 
phase  adhere  to  the  surface  of  solids,  was  not  included  in  the  original 
model.  Since  the  net  effect  of  this  process  would  be  to  remove  nitrogen 
and  phosphorus  from  the  water  column,  where  predictions  were  high,  the 
process  was  programmed  and  evaluated  in  CE-QUAL-R1.  The  first  algorithm 
had  a  constant  amount  of  phosphorus  and  nitrogen  adsorbed  per  gram  of 
suspended  solids,  which  included  detritus,  algae,  and  suspended  solids. 
The  phosphorus  and  nitrogen  were  moved  to  the  sediment  or  to  next  lower 
layer,  depending  on  the  average  settling  rate  of  the  solids. 
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where 


A  =  amount  of  nutrient  settled  from  a  layer 
due  to  adsorption  (g/layer/tiroe  step) 

3 

V  =  volume  of  water  in  the  layer  (m  /layer) 

K  =  coefficient  representing  the  amount  of 
nutrient  adsorbed  per  unit  of  suspended 

solids  ( 8 

\  g  solids/m  / 

C  Cd.  C  j,  and  C ^  =  the  concentration  of  suspended  solids, 

detritus,  and  the  two  algal  compartments 
(g/m  ) 

S  ,  S  ,  S  .,  and  S  =  the  settling  rate  of  the  above  compart- 
S  3  3  ments  (m/time  step) 

T  =  the  layer  thickness  (m) 

34.  There  were  problems  with  this  formulation  because  it  allowed 
most  of  the  inorganic  nutrients  to  be  settled  in  layers  with  low  nutri¬ 
ent  concentrations,  while  at  the  same  time  not  having  the  desired  effect 
on  layers  with  high  nutrient  concentrations.  To  resolve  this  problem, 
the  coefficient  K&  in  Equation  1  was  made  a  function  of  the  nutrient 
concentration  (after  Hwang,  Lackie,  and  Huang  1976). 
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where 
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N  =  inorganic  nutrient  concentration  (g/m  ) 

K  =  maximum  amount  of  solute  per  amount  of  solids  (g  nutrient/g 
m  solids) 

3 

=  adsorption  coefficient  fl/(time  step  g/m  )]  divided  by  the 
desorption  coefficient  (1/time  step) 

35.  Variable  was  read  in  as  a  single  coefficient  since  equi¬ 

librium  concentrations  were  assumed.  Equation  2  is  known  as  the  Lang¬ 
muir  isotherm,  which  has  been  used  to  describe  nutrient  adsorption  (Ku, 
Di  Giano,  and  Feng  1978) .  This  formulation  noticeably  improved  results 
and  was  therefore  considered  satisfactory. 


Addition  of  Refractory  Dissolved  Organics 

36.  The  original  model  included  one  variable  for  dissolved 
organic  matter.  Input  to  this  compartment  was  from  upstream  flow,  while 
output  included  decomposition  and  flow  out  of  the  system.  Maximum  decay 
was  specified  by  a  coefficient  and  reduced  as  a  function  of  temperature. 
This  representation  for  decomposition  of  dissolved  organic  matter  was 
too  simple.  The  dissolved  organic  matter  compartment  in  a  reservoir 
actually  consists  of  a  number  of  compounds  that  decompose  at  varying 
rates.  Rates  near  0.4/day  have  been  measured  for  labile  compounds,  and 
near  0.004/day  for  refractory  compounds. 

37.  Many  of  the  labile  compounds  in  a  reservoir  come  from  algae 
through  processes  of  mortality  and  respiration.  Materials  coming  from 
upstream  usually  have  been  undergoing  decomposition  in  the  river  and  are 
more  refractory.  To  better  represent  these  processes,  two  dissolved 
organic  matter  compartments  were  included  in  the  model:  one  represent¬ 
ing  refractory  organic  compounds  and  the  other  representing  labile  com¬ 
pounds.  Each  of  the  compartments  has  a  different  decomposition  rate. 

As  the  labile  compartment  decays,  some  of  the  products  go  to  the  inor¬ 
ganic  nutrients  and  some  to  the  refractory  organic  compartment.  A  dia¬ 
gram  of  these  compartments  is  included  as  Figure  6.  These  changes 
improved  model  predictions. 

Inflow  Modifications 


38.  CE-QUAL-R1,  being  one-dimensional,  does  not  allow  the  model 
to  predict  variations  as  one  moves  from  the  headwater  to  the  damsite  of 
a  reservoir.  This  could  cause  problems,  because  the  concentrations  of 
inflowing  constituents  are  usually  measured  at  or  above  the  headwater 
area,  and  the  model's  predictions  are  best  suited  near  the  deepest  part 
of  the  pool,  which  usually  occurs  near  the  dam.  Any  processes  occurring 
between  the  headwater  and  dam  which  can  change  the  concentrations  of  in¬ 
flowing  constituents  would  not  be  represented  and  can  lead  to  erroneous 
results.  Since  phosphorus  and  nitrogen  can  be  taken  up  by  phytoplankton 


or  adsorbed  to  particulate  materials  and  settled  to  the  headwater  bot¬ 
tom,  the  concentration  of  inorganic  nutrients  entering  the  main  part  of 
the  reservoir  may  have  been  lower  than  used  for  the  DeGray  simulations. 
To  demonstrate  the  effect  this  might  have  had  on  the  high  concentration 
of  phosphorus  and  nitrogen  in  the  hypolimnion  as  predicted  by  the  model, 
the  concentration  of  nutrients  in  the  inflow  was  lowered  by  70  percent 
and  100  percent.  This  change  had  virtually  no  effect  on  the  high  hypo- 
limnetic  concentrations  of  phosphorus  and  nitrogen,  but  did  have  an 
adverse  effect  on  the  prediction  of  algae  in  the  spring  in  the  epilim- 
nion.  For  this  reason  no  permanent  change  was  made  to  the  model  con¬ 
cerning  inflowing  concentrations  of  nutrients. 

Phytoplankton  Modifications 


39.  The  Monod  equation  for  photosynthesis  as  a  function  of  light 
was  replaced  with  Steele's  (1962)  light  equation. 
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where 

F  =  fraction  of  maximum  photosynthesis 

2 

I  =  average  light  for  a  particular  layer  (kcal/m  /hr) 

K  =  Steele's  light  saturation  coefficient  or  the  amount  of 
light  (kcal/m2/hr)  at  maximum  photosynthesis 

Steele's  equation  is  more  widely  accepted  and  can  represent  photoinhibi¬ 
tion  at  irradiance  levels  greater  than  the  light  saturation  (Kg)  level. 

40.  Loss  terms  for  photorespiration,  or  excretion,  and  mortality 
were  also  added  to  the  phytoplankton  compartment.  Photorespiration  re¬ 
sults  from  the  oxidation  of  ribulose  diphosphate  instead  of  its  carboxy- 
lation  to  yield  glycolate.  This  process  has  been  shown  to  be  most  sen¬ 
sitive  to  low  and  high  irradiance  levels  and  can  be  represented  by  the 
function 
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E  =  fraction  of  the  compartment  lost  to  excretion  (per  time 
step) 

=  maximum  excretion  rate  (per  time  step) 
and  F  is  calculated  in  Equation  3. 

Photorespiration  rates  increase  at  both  very  high  and  very  low  intensi¬ 
ties.  Extracellular  release  of  photoassimilated  carbon  can  range  from  a 
baseline  of  2  percent  to  as  high  as  40  percent. 

41.  Another  change  dealt  with  the  products  of  algal  respiration. 
Based  on  constant  stoichiometry  of  algae,  products  from  total  respira¬ 
tion  were  fractioned  into  the  compartments  of  phosphorus,  inorganic  car¬ 
bon,  and  ammonia.  Thus,  the  respired  biomass  was  again  immediately 
available  for  use  by  the  algae.  In  actuality,  this  cycling  scheme  is 
erroneous  and  can  create  problems  for  phosphorus  and  ammonia  dynamics. 
The  major  products  of  respiration  are  carbon  dioxide  and  organic  com¬ 
pounds.  These  organic  compounds  are  later  remineralized  into  readily 
available  forms  of  nitrogen,  phosphorus,  and  inorganic  carbon.  To  rep¬ 
resent  this  process  correctly,  the  products  of  algal  respiration  should 
be  carbon  dioxide  and  labile  dissolved  organics.  This  would  alleviate 
the  problem  of  immediate  recycling,  since  dissolved  organics  decay  at  a 
slower  rate.  Unfortunately,  this  can  cause  a  mass  balance  problem  for 
carbon  because  of  the  model  simplification  of  constant  stoichiometry. 

Tn  this  case,  the  same  carbon  would  be  released  twice,  once  as  respired 
carbon  dioxide  and  once  when  dissolved  organic  matter  decomposes.  In 
order  to  maintain  a  continuity  of  mass  and  to  slow  the  immediate  recycl¬ 
ing  of  nutrients,  the  products  of  photorespiration  increase  the  labile 
dissolved  organic  matter  compartment,  while  those  from  dark  respiration 
will  continue  to  increase  the  NH4-N,  inorganic  carbon,  and  P04-P 
compartments . 

42.  The  mortality  term  added  to  the  algae  compartment  represents 
a  loss  rate  due  to  the  senescence  of  algae.  It  is  a  difficult  parameter 
to  measure  in  the  field,  but  is  generally  regarded  as  being  less  than 


10  percent  of  the  total  loss  terms.  This  rate  is  accelerated  by 
critically  high  temperatures  and  can  be  simulated  by  the  equation 
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where 

D  =  fraction  of  the  compartment  that  dies  (per  time  step) 

K  =  maximum  mortality  fraction  (per  time  step) 

T  =  layer  temperature  (°C) 

T  =  temperature  at  which  total  mortality  occurs  (°C) 
m 

Mortality  is  fractioned  between  the  detritus  and  labile  dissolved  or¬ 
ganic  matter  compartments . 

A3.  In  the  past,  no  photosynthesis  was  allowed  to  occur  below  the 
1-percent  light  level.  This  is  a  somewhat  arbitrary  standard  that  indi¬ 
cates  the  depth  where  only  1  percent  of  the  incident  irradiance  is  pene¬ 
trating  from  the  surface.  This  has  been  referred  to  as  the  compensation 
depth,  a  point  where  phytoplankton  respiration  exceeds  photosynthesis. 
Photosynthesis  does  not  stop  here,  however.  Futhermore,  shade-adapted 
species  are  known  to  acclimate  to  low  irradiance  levels  by  reducing 
their  respiration  rates,  thereby  increasing  the  photosynthesis/respira¬ 
tion  ratio  compensation  level.  This  can  be  very  important  in  some  sys¬ 
tems  and  under  conditions  of  stratification.  The  elimination  of  the 
condition  has  been  tested  in  the  model  under  the  assumption  that  respi¬ 
ration  processes  will  balance  out  with  reduced  photosynthetic  ability. 
The  combined  changes  to  the  algal  subroutine  produced  improved  predic¬ 
tions.  All  of  the  structural  changes  to  the  phytoplankton  compartment 
are  shown  in  Figure  7. 
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Fluxes 


44.  Initial  calibration  of  CE-QUAL-Rl  for  DeGray  in  1979  started 
with  coefficients,  estimated  by  Johnson  and  Ford  (1981),  for  the  thermal 
model  CE-THF.RM-R1  and  textbook  values  for  biological  and  chemical 
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coefficients.  The  initial  simulation  showed  decreases  in  the  oxygen 
concentrations  in  the  metalimnion  and  the  hypolimnion,  although  at  a 
faster  rate  than  observations  indicated.  Instead  of  arbitrarily  chang¬ 
ing  coefficients  until  satisfactory  concentrations  were  predicted,  the 
fluxes  predicted  by  the  model  were  compared  with  field  measurements  or 
literature  values.  Although  flux  information  from  the  model  was  avail¬ 
able  on  a  per  layer  basis,  the  study  concentrated  on  fluxes  at  the  res¬ 
ervoir  level  to  search  for  gross  discrepancies. 

45.  At  the  time  the  flux  information  was  being  used  to  help  cali¬ 
brate  the  1979  data  set,  there  were  3  processes,  on  a  reservoir  basis, 
that  could  increase  oxygen  concentrations  and  13  processes  that  could 
decrease  oxygen  concentrations  (Table  1).  Of  the  13  that  resulted  in 
losses,  3  processes  accounted  for  80  percent  of  the  oxygen  loss.  The 
model  predicted  that  40  percent  of  the  oxygen  loss  was  used  in  algal 
respiration,  29  percent  was  lost  in  the  decomposition  of  dissolved 
organic  matter,  and  11  percent  was  due  to  the  decomposition  of  sedi¬ 
ments.  In  addition,  algal  respiration  accounted  for  62  percent  of  pri¬ 
mary  production.  Carbon  primary  production  data  were  available  for 
heGray  in  1979.*  The  samples  were  collected  once  per  month  and  reported 
as  milligrams  of  carbon  per  square  meter  per  hour.  Assuming  constant 
values  for  the  day  and  month,  primary  productivity  would  be  approxi- 
matelv  4  x  10  kg  carbon  for  the  reservoir  for  the  simulation  period. 

The  model  value  for  the  period  was  1.5  x  10^  kg  of  carbon,  which  is 
approximately  3.7  times  the  measured  value.  In  addition,  Whittaker's 
(1°75)  general  figure  for  algal  respiration  of  30  to  40  percent  of  gross 
primary  production  indicates  the  model  prediction  of  62  percent  was  high. 

Ub.  No  data  were  available  for  DeGray  Take  concerning  decomposi¬ 
tion  of  dissolved  organic  matter  or  sediments.  Tn  studies  on  other  res- 

i 

ervoirs,  58  to  176  mg  oxvgen/m  /day  were  utilized  bv  the  sediments  (Gun¬ 
nison,  Chen,  and  Brannon  1983).  The  initial  value  predicted  by  the 

a 

model  was  475  mg  oxygen/mVday .  Since  the  major  process  concerning 
dissolved  organic  matter  was  decomposition  and  since  predicted 
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*  Personal  Communication,  R.  H.  Kennedy,  op.  cit. 
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concentrations  were  low,  the  decomposition  coefficient  was  probably  high 
and  was  therefore  lowered.  Although  collected  in  1980,  sediment  trap 
data  were  also  available  for  DeGray  Lake.  This  information  gave  an 
indication  of  the  amount  of  algae  settling  at  three  depths:  5,  15,  and 
45  m.*  Data  were  collected  at  approximately  monthly  intervals  for  six 
periods  between  February  and  August  and  were  reported  as  milligrams 
chlorophyll  a  settling  per  square  meter  per  day.  Assuming  the  average 
rate  for  that  period  remained  constant  for  the  entire  simulation,  and 

3 

using  Spangler's  (1969)  conversion  factor  of  0.23  g/m  dry  weight  equals 
milligrams  of  chlorophyll  a  per  cubic  meter,  the  algal  mass  settling  at 
these  depths  was  calculated.  The  figures  were  3.8  x  10^,  3.4  x  10^,  and 

4 

1.0  x  10  kg/reservoir/339  days  for  the  depths  of  5,  15,  and  45  m, 

respectively.  Corresponding  predictions  were  4.2  x  10^,  9.8  x  10"*,  and 
4 

1.8  x  10  ,  respectively. 

47.  With  this  evidence  in  mind,  coefficients  dealing  with  algal 
production,  respiration  and  settling,  and  sediment  and  dissolved  organic 
matter  decomposition  were  reduced  in  order  to  bring  model  predictions  in 
line  with  measured  flux  values.  During  the  ensuing  calibration  simula¬ 
tions,  as  the  fluxes  were  brought  more  in  line  with  measured  values,  the 
oxygen  concentrations  also  improved.  The  predicted  flux  values  for  the 
beginning  and  ending  of  the  calibration  exercise  are  listed  in  Table  6. 
Predicted  fluxes  for  primary  production,  algal  respiration,  and  sediment 
oxygen  demand  were  near  measured  values.  For  algal  settling  at  5  and 
15  m,  predictions  were  not  as  good  as  the  other  fluxes,  but  did  show  im¬ 
provement  over  the  initial  simulations.  Since  predicted  settling  at  5  m 
was  higher  and  at  15  m  was  lower  than  measured  values,  predictions  were 
considered  satisfactory.  The  predicted  value  for  algal  settling  at  45  m 
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PART  V:  CONCLUSIONS 


48.  CE-QUAL-R1  was  evaluated  with  data  collected  at  DeGray  Lake, 
Ark.  Data  collected  in  :979  were  used  for  calibration,  while  data  col¬ 
lected  in  1980  were  used  for  confirmation.  Evaluation  consisted  of 
tests  of  the  code  and  comparison  of  measured  values  to  model  predic¬ 
tions.  Tests  to  evaluate  the  software  pinpointed  a  number  of  errors  in 
the  original  model.  Corrections  ranged  from  changing  misspelled  vari¬ 
able  names  to  replacing  algorithms  or  changing  part  of  the  solution 
scheme . 

49.  Evaluation  of  output  from  the  original  model  showed  that  pre¬ 
dicted  ortho-P  and  ammonia-nitrogen  values  were  not  at  all  satisfactory, 
especially  in  the  hvpolimnion  during  the  latter  months  of  the  year.  The 
problem  persisted  through  all  calibration  and  confirmation  simulations, 
indicating  that  the  model  structure  may  not  have  been  correct.  A  number 
of  changes  to  algorithms,  processes,  and  compartments  were  programmed 
and  tested  in  order  to  simplify  the  model  or  improve  predictions. 

Changes  included  uniting  some  compartments,  adding  a  dissolved  refrac¬ 
tory  organic  compartment,  adding  the  process  of  adsorption,  and  changing 
algorithms  dealing  with  phytoplankton.  The  model  has  been  improved  in 
that  it  requires  22  fewer  coefficients,  while  predictions  more  closely 
match  measured  values,  especially  for  ortho-P  and  ammonia-nitrogen.  For 
the  confirmation  simulations  for  which  3,536  comparisons  were  made,  the 
average  value  for  the  Reliability  Index  improved  from  3.49  to  2.59. 
Graphs  of  model  output  before  (original  model)  and  after  revisions 
(final  model)  are  shown  in  Figures  1  and  2. 

50.  Model  changes  based  on  the  results  of  this  and  other  studies 
are  being  incorporated  in  a  new  version  of  the  CE-QUAL-R1  User's  Manual. 
Coefficients  for  the  final  simulations  are  listed  in  Appendix  A;  Appen¬ 
dix  B  presents  the  initial  values  of  state  variables.  Complete  data 
sets  for  1979  and  1980  can  be  found  in  Appendixes  C  and  D,  respectively. 
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Table  1 

Process  Interactions  for  Variables  in  CE-OUAL-R J ,  Original 


Version  (Environmental  Laboratory  1982a) 
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Definitions  of  abbreviations  are  as  follows: 

D  -  Decay  or  decomposition 
DOM  ■  Dissolved  organic  matter 
E  -  Egestlon 

F  •  Diffusion  and  convection 
G  ■  Gain  or  loss  caused  by  layer  depth  change 
H  *  Fishing  harvest 
I  -  Ingestion 

N  -  Nonpredatorv  mortality 
0  -  Outflow 
P  •  Photosynthesis 
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R  ■  Respiration 
S  -  Settling 
SS  *  Suspended  solids 
TDS  ■  Total  dissolved  solids 
W  -  Inflow 

X  -  Exchange  at  the  air-water 
interface 

Y  ■  Settling,  diffusion  and  convection 
Z  -  Egestion  and  nonpredatory 
mortality 
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Statistical  Results  from  the  Confirmation  Simulation  of  DeGray  Lake,  1980 
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Table  6 

arlsons  of  Initial  and  Final  Simulations  to  Measured  Values 
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Figure  1.  Model  predictions  versus  measured  values  (cir¬ 
cles)  for  DeGray  Lake  in  1979.  The  dashed  line  represents 
predictions  from  the  original  model;  the  solid  line  repre¬ 
sents  those  from  the  final  mode]  (Sheet  1  of  4) 
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Figure  1.  (Sheet  3  of  4) 
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Figure  2.  Model  predictions  versus  measured  values  (cir¬ 
cles)  for  DeGray  Lake  in  1980.  The  dashed  line  represents 
predictions  from  the  original  model;  the  solid  line  repre¬ 
sents  those  from  the  final  model  (Sheet  1  of  4) 
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Figure  2.  (Sheet  3  of  4) 
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Figure  3.  Interaction  of  fish  with  other  compartments 
after  model  improvements 
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Figure  4.  Interaction  of  benthos  plus  sediment  with  other 
compartments  after  model  improvements 


H 


OXYGEN 


INFLOW 


- 1 -  PHOTOSYNTHESIS 

NITRITE  +  NITRATE  - ► 

NITROGEN  OUTFLOW 


EACH  ALGAL 
GROUP 


|  NH4-N 


DECAY 


1 - 1 - 3 

' 

- * 

DIFFUSION  | 

OTHER  N02-N  + 

N03-N  LAYERS 

Figure  5.  Interaction  of  nitrite  plus  nitrate-nitrogen  with  other 
compartments  after  model  improvements 


Figure  6.  Interaction  of  the  two  dissolved  organic  matter  compartments 
with  other  variables  after  model  improvements 


Figure  7.  Interaction  of  each  algal  group  with  other 
compartments  after  model  improvements 


_ Parameter 

Reservoir  description 


Model 

icronvn 


Units 


Value 


Number  of  outlets 

NOUTS 

3 

Number  of  tributaries 

NTRIBS 

2 

-  -  j 

Latitude 

XLAT 

34.2 

-*->/ 

Longitude 

XLON 

decimal  degrees 

93.1 

•  • 

v';/, 

Reservoir  length 

RLEN 

decimal  degrees 

13000.0 

Minimum  layer  thickness 

SDZMIN 

m 

0.5 

Maximum  layer  thickness 

SDZMAX 

m 

2.0 

Area  coefficient 

ACOEF(l) 

m 

561.81 

Area  coefficient 

AC0EF(2) 

2.79 

Width  coefficient 

WCOEF(l) 

47.7 

rjk 

Width  coef  f ic ient 

WC0EF(2) 

0.55 

y 

Port  1  elevation 

ELOUT(l) 

m 

56.4 

*,.**,•» 

Port  1  area 

AROUT ( 1 ) 

2 

m 

31.2 

!yv* 

Port  2  elevation 

EL0UT(2) 

m 

51.8 

.«* 

Port  2  area 

AROUT (2) 

2 

m 

31.2 

Port  3  elevation 

EL0UT(3) 

m 

44.4 

Port  3  area 

AROUT (3) 

2 

m 

31.2 

•N-S- 

Physical  coefficients 
Turbidity  factor 
Wind  coefficient 
Wind  coefficient 
Sheltering  coefficient 


TURB 

AA 

BB 

SHELCF 


m/ (sec 
1/mb 


Penetrative  convection  fraction  PEFRAC 


Wind  mixing  coefficient 
Advection  mixing  coefficient 
Critical  density 


CDIFW 

CDIFF 

CDENS 

(Continued) 


2.0 

0 

0.12  x  10-8 
1.0 
0.01 
0.004 
0.0007 
2.0 
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Extinction  coefficient 


For  water 

EXCO 

1  /m 

0.45 

For  inorganic  solids 

EXTINS 

1/m 

0.01 

For  organic  solids 

EXTINP 

1/m 

0.1 

Surface  radiation  fraction 

SURFRAC 

0.4 

Reaeration  coefficient  -  oxygen 

DM02 

2, 

m  /sec 

2.04  x  10“ 

Reaeration  coefficient  -  CO^ 

DMC02 

2, 

m  /sec 

1.63  x  10" 

Stoichiometry 

Carbon  -  organics 

0.46 

Nitrogen  -  organics 

0.05 

Phosphorus  -  organics 

0.004 

Oxygen  -  ammonia 

02NH3 

4.57 

Oxygen  -  nitrite  +  nitrate 

02N02 

1.14 

Oxygen  -  detritus  decay 

02DET 

1.4 

Oxygen  -  respiration 

02RESP 

1.1 

Oxygen  -  photosynthesis 

02  FAC 

1.4 

Oxygen  -  dissolved  organic  decay 

02D0M 

1.4 

Oxygen  -  manganese 

02MN2 

0.15 

Oxygen  -  iron 

02FE2 

0.14 

Oxygen  -  sulfide 

02S2 

2.0 

Organic 

Phytoplankton 

Respiration  rate 

TPRESP 

1/day 

0.017 

Gross  production  rate 

TPMAX(l) 

1/day 

1.1 

Fraction  death  to  detritus 

ALDIGO 

0.25 

Settling  rate 

TSETL(l) 

m/day 

0.14 

(Continued) 


Parameter 

Model 

Acronym 

Units 

Value 

Organic  (Cont.) 

Half-saturation 

Carbon 

PS2C02(1) 

mg/£ 

0.12 

Nitrogen 

PS2N( 1) 

mg  /*• 

0.014 

Phosphorus 

PS2P04 (1) 

mg/ £ 

0.009 

Light  saturation  level 

PISAT(l) 

kcal/m^/hr 

50.0 

Excretion  rate 

TPEXCR(l) 

1/day 

0.01 

Mortality  rate 

TPMORT ( 1 ) 

1/day 

0.01 

Temperature  multipliers 

Low  threshold 

ALG1T1 

°C 

0 

Low  optimum 

ALG1T2 

°C 

26.0 

High  optimum 

ALG1T3 

°C 

30.0 

High  threshold 

ALG1T4 

°C 

35.0 

Low  minimum 

ALG1K1 

0.1 

High  minimum 

ALG1K4 

0.1 

Zooplankton 

Ingestion  rate 

TZMAX 

I /day 

0.44 

Mortality  rate 

TZMORT 

1/day 

0.01 

Efficiency 

ZEFFIC 

0.5 

Food  preference 

For  algae  1 

PREF(l) 

0.5 

For  algae  2 

PREF(2) 

0 

For  detritus 

PREF(3) 

0.5 

Respiration  rate 

TZRESP 

0. 14 

Half-saturation 

ZS2P 

mg/ £ 

0.3 

Temperature  multiplier 

Low  threshold 

Z00T1 

°C 

0 

Low  optimum 

Z00T2 

°C 

20.0 

(Continued) 


Parameter 

Model 

Acronym 

Units 

Value 

Organic  (Cont.) 

High  optimum 

Z00T3 

°C 

26.0 

High  threshold 

Z00T4 

°C 

36.0 

Low  minimum 

Z00K1 

0.1 

High  minimum 

Z00K4 

0.1 

Fish 

Ingestion  rate 

TFMAX 

1/day 

0.015 

Half-saturation 

FS2FSH 

0.2 

Food  preference 

For  benthos  and  sediment 

FPSED 

0.03 

For  algae  1 

FPALG(l) 

0.37 

For  algae  2 

FPALG(2) 

0 

For  zooplankton 

FPZOO 

0.34 

For  detritus 

FPDET 

0.26 

Efficiency 

FEFFIC 

0.8 

Mortality  rate 

TFMORT 

1/day 

0.01 

Respiration  rate 

TFRESP 

1/day 

0.01 

Temperature  multiplier 

Low  threshold 

FSHIT1 

°C 

1.0 

Low  optimum 

FSH1T2 

°C 

24.4 

High  optimum 

FSH1T3 

°C 

28.4 

High  threshold 

FSH1T4 

°c 

35.2 

Low  minimum 

FSH1K1 

0.1 

High  minimum 

FSH1K4 

0.1 

Decomposition 

Labile  organics 

TDOMDK 

1/day 

0.032 

Ammonia 

TNH3DK 

1  /  day 

0.08 

(Continued) 


Parameter 


Organic  (Cont.) 


Model 

Acron 


Units 


Value 


•  V  V  V  ■ 


Detritus 

TDETDK 

1/day 

0.009 

-V-V-V-W-*- 

Coliforms 

TCOLDK 

1/day 

1.4 

'.V- V.\V>V. 

Sediment 

TSEDDK 

1/day 

0.008 

Refractory  organics 

TRFRDK 

1/day 

0.005 

Labile  to  refractory 

TDOMRF 

1  /day 

0.005 

Temperature  multipliers 

DOM  low  threshold 

D0MT1 

°C 

2.0 

DOM  optimum 

D0MT2 

°C 

20.0 

DOM  low  minimum 

D0MK1 

0.12 

NH3  low  threshold 

NH3T1 

°C 

2.0 

. v- 

NH3  optimum 

NH3T2 

°C 

32.0 

NH3  low  minimum 

NH3K1 

0.1 

mi-M 

N02  low  threshold 

N02TL 

6C 

2.0 

NO 2  optimum 

N02T2 

°C 

32.0 

•  ■  «*■.  •*. 

NO 2  low  minimum 

N02K1 

0.1 

f  1 

Inorganics 

Solids  settling 

TSSETL 

m/day 

0.05 

P04  adsorption 

ADSRBP 

1/m3 

150.0 

Nitrogen  adsorption 

ADSRBN 

1/m3 

125.0 

P04  adsorption 

ADMAXP 

g/g 

0.007 

Nitrogen  maximum  adsorption 

ADMAXN 

g/g 

0.005 

APPENDIX  B:  INITIAL  VALUES  OF  STATE  VARIABLES  FOR  THE 
1979  AND  1980  DEGRAY  SIMULATIONS 


1979  1980 


Variables 

Units 

Lowest 

Highest 

Lowest 

Highest 

Temperature 

°C 

5.2 

5.5 

6.4 

9.0 

Oxygen 

mg/ 5 

9.2 

9.5 

0.8 

9.9 

Algae 

mg/ 5 

0.0 

1.6 

0.0 

1.0 

Zooplankton 

mg/ 5 

0.0 

0.01 

0.001 

0.001 

Coliforms 

100/5 

20.0 

1060.0 

0.0 

2.0 

Ammonia  -  N 

mg/  5 

0.0 

0.008 

0.01 

0.06 

N02-N  +  N03-N 

mg/ 5 

0.09 

0.23 

0.06 

0.33 

P04-P 

mg/  J l 

0.001 

0.004 

0.0 

0.001 

Detritus 

mg  /  i 

0.0 

0.0 

0.2 

0.7 

Sediment 

g/m2 

101.0 

101.0 

100.0 

100.0 

Alkalinity 

mg/  it 

14.8 

20.0 

15.0 

115.0 

Total  dissolved  solids 

mg/ 5 

49.0 

69.0 

0.0 

32.0 

Suspended  solids 

mg/5 

6.0 

6.0 

6.0 

28.0 

Labile  organics 

mg/ 5 

1.9 

4.4 

1.7 

3.1 

Refractory  organics 

mg/5. 

7.8 

17.6 

5.8 

12.5 

Inorganic  carbon 

mg/ 5. 

4.0 

8.7 

7.0 

53.0 

Carbon  dioxide 

mg/ 5 

0.46 

3.9 

3.4 

26.0 

pH 

mg/  5 

6.6 

7.4 

5.9 

6.5 

Particulate  manganese 

mg/ 5. 

0.1 

0. 1 

0.0 

0.0 

Sediment  manganese 

mg/ 5 

30.0 

30.0 

30.0 

30.0 

Dissolved  manganese 

mg/ £ 

0,0 

0.0 

0.  1 

1.2 

Particulate  iron 

mg/ 5, 

0.0 

0.  I 

0.0 

0.0 

Sediment  iron 

mg/ 5. 

600.0 

600.0 

600.0 

600.0 

Dissolved  iron 

mg  / 1 

0.0 

0.0 

0.  1 

0.7 

Iron  sulfide  -  sediment 

mg/ 5 

0.0 

0.0 

0.0 

0.0 

Iron  sulfide  -  water 

mg/5 

0.0 

0.0 

0.0 

0.0 

Sulfate 

mg/ 5. 

3.0 

4.0 

3.0 

4.0 

Sediment  sulfur 

mg/5. 

10.0 

10.0 

10.0 

10.0 

Sulfide 

mg/5 

0.0 

0.0 

0.0 

0.  1 

Sediment  P 

mg/ 5 

20.0 

20.0 

20.0 

20.0 

Sediment  N 

mg/  5 

50.0 

50.0 

50.0 

50.0 

APPENDIX  C:  1979  DATA  SET  FOR  DEGRAY  LAKE 


TITLE 

T'TLE 


DEGRAY  1979 
FINAL  CALIBRATION 


TITLE 

TITLE 

TITLE 


JOB 

1 

362 

24 

720 

25 

79 

1 

OUTPUT 

COMPLETE 

PHYS1 

3 

2 

60 

34.2 

93.1 

2 

0 

PHYS2 

13000 

.  5 

2.0 

PHYS2+ 

1.25 

1.25 

1  . 

1  . 

1  . 

1 . 

1 . 

PHYS2+ 

1  . 

1  . 

1 . 

1 . 

1  . 

1  . 

1 . 

PHYS2+ 

1  . 

1  . 

1. 

1  . 

1  . 

1  . 

1 . 

PHYS2* 

1  . 

1  . 

1  . 

1 . 

1 . 

1 . 

1 . 

PHYS2+ 

1  . 

1  . 

1  . 

1  . 

1  . 

1  . 

1 . 

PHYS2+ 

1  . 

1  . 

1 . 

1 . 

1  . 

1  . 

1 . 

PHYS2+ 

1  . 

1 . 

1 . 

1  . 

1 . 

1  . 

STRUCT 

PORT 

CHOICE 

SPECIFIED 

PHYS3 

56.4 

5.58 

5.58 

PHYS3 

51.8 

5.58 

5.58 

PHYS3 

44.4 

5.58 

5.58 

PHYS4 

561.81 

2.79 

PHYS5 

47.70 

0.55 

MIXING 

1.0 

.01 

.00004 

.000008 

2.0 

LIGHT 

.45 

0.4 

.01 

DIFC2 

2.04-09  1. 

.63-09 

ALG1 

.10 

0.017 

.25 

ALG2 

1.10 

0.14 

.009 

.014 

0.12 

50  . 

.010 

ALG3 

0.8 

0.14 

.009 

.01 

0  .  1 

20  . 

.  020 

ALG3A 

1.0 

0 .14 

.009 

.01 

0.1 

54. 

.001 

ALG3++ 

.05 

ALGA 

0 

26 

30 

35 

0.1 

0.1 

ALG5 

4 

26 

36 

40 

0  .  1 

0.1 

ALG5  + 

2 

26 

32 

37 

0.1 

0.1 

PL  AN  T 1 

1.2 

.2 

.  1 

.05 

.4 

.3 

.  3 

PL  ANT  2 

.02 

.05 

.01 

.005 

10. 

30  . 

.5 

PLANTS 

2. 

25. 

29. 

38. 

.  1 

.  1 

ZOOl 

.44 

.010 

0 .50 

0.5 

0.0 

0.0 

0.5 

Z002 

.30 

0  .  0 

20 

26 

36 

0.1 

0  .  1 

DET1 

.35 

4.0 

28 

0.01 

FISH1 

.0150 

.2 

.03 

.  37 

0  . 

0.0 

.34 

FISH2 

1  . 

24.4 

28.4 

35.2 

.  1 

.  1 

.8 

DECAY1 

0.032 

0  .  08 

0.009 

1.4 

.  0080 

.  0050 

.005 

DECAY2 

2 

20 

.12 

DECAYS 

2 

32 

0.1 

DECAY4 

2 

32 

0.1 

SSETL 

.05 

150. 

125. 

.007 

.005 

TMP 

1 . 04 

CHEM 

4 . 57 

1 .14 

1.4 

1 .  1 

1.4 

1.4 

0.15 

ANAER1 

0.5 

5.0 

ANAER2 

0.05 

0.02 

0 

5 

35 

40 

0.1 

ANAER3 

0.10 

0 

5 

35 

40 

0  .  1 

0  .  1 

ANAER4 

0.00 

0 

5 

35 

40 

0  .  1 

0  .  1 

ANAER5 

0 .05 

0.03 

0 

5 

35 

40 

0  .  1 

ANAER6 

0.10 

0 

5 

35 

40 

0  .  1 

0  .  1 

ANAER7 

0.00 

0 . 0 

0 

5 

35 

40 

0  .  1 

lANAERB 

0 .90 

0 

5 

35 

40 

0.1 

0  .  1 

ANAER9 

0.50 

0  .  5 

0 

5 

35 

40 

0  .  1 

ANAER10 

0.001 

0 

•5 

35 

40 

0 . 1 

0  .  1 

ANAER11 

0.000012 

0 

5 

35 

40 

0  .  1 

0  . 1 

s  \ 


.001 

1 . 

.0001 

.22 

5.3 

60.4 

)  0.0 

4.0 

.001 

1  . 

.000  1 

.21 

5 . 2 

59.6 

1  0.0 

4 . 0 

.001 

1  . 

.000  1 

.  21 

5.3 

58.8 

1  0.0 

4.0 

.001 

1  . 

.  0  0  0  1 

.21 

5.3 

58. 

1  0.0 

4  .  0 

.001 

1  . 

.000  1 

.21 

5.  3 

60.7 

t  0.0 

4 . 0 

n  n  i 

i 

Trn 


■v  _ 

%  ^ 

"F  *  "  -  ' 

INIT5 

12000  . 

0.0 

200  . 

400  . 

1000  . 

.001 

INIT2 

57.5 

1.42 

0. 

17. 

1.001 

.  0001 

INIT3 

9.9 

9.4 

5.4 

INIT4 

6  . 

mmm 

0.0 

INIT5 

12000. 

BH 

■ETTTjB 

.  001 

INIT2 

58.5 

.8 

Wmm 

16.5 

1.001 

.005 

INIT3 

.0002 

11 . 

9.4 

101.1 

5.4 

INIT4 

6  . 

0  . 1 

0  .  0 

INIT5 

12000  . 

0.0 

200  . 

400  . 

1  00  0  . 

.001 

I  HI  12 

59.5 

1  .  5 

0. 

15.  7 

1  .  001 

.  007 

INIT3 

.  0002 

13.7 

9 . 5 

.002 

101.1 

5 . 4 

I N I  T  4 

6  . 

0.1 

0  .  0 

0.0 

0  .  0 

0  .  0 

INIT5 

12000 . 

0  .  0 

200  . 

400  . 

100  0. 

.001 

INIT2 

59.7 

1.4 

0  . 

14 .8 

1.001 

.  008 

INIT3 

.0002 

16  .  3 

9.5 

.002 

101.1 

5.5 

IHIT4 

6  . 

0.1 

0  .  0 

0  .  0 

0  .  0 

0  .  0 

i  h  i  r  5 

12000  . 

0.0 

200  . 

400  . 

1000  . 

.001 

INIT  PL 

0  . 

0  . 

0  . 

0  . 

0  . 

0. 

FILES  PLTWC  PLDG15  PLDG16  PLDG17  FLUX 
ID  DEGRAY7 9  MOV  2  83  HI  CDC  +  JU 


WEATH1 

24 

365 

W2 

LTLRC 

79 

1  1 

0  .  9 

-2.0 

-7.4 

1013.5 

19.7 

W2 

LTLRC 

79 

1  2 

0  .  0 

-7  .2 

-17.2 

1022.9 

16.2 

U2 

LTLRC 

79 

1  3 

0.3 

-3.3 

-11.9 

1027.1 

5.8 

Ul  2 

LTLRC 

79 

1  4 

0.6 

-0  .  3 

-7  .8 

1028.3 

6  .  3 

W2 

LTLRC 

79 

1  5 

1  .  0 

0 . 3 

-3.5 

1  023.1 

10  .4 

l**2 

LTLRC 

79 

1  6 

1  .  0 

-2.6 

-6.5 

1015.1 

13.7 

1*12 

LTLRC 

79 

1  7 

0  .  9 

-1.5 

-6  .  3 

1011.5 

13.0 

1*12 

LTLRC 

79 

1  8 

0  .  1 

-6  .  3 

-13.2 

1022.1 

8. 1 

W2 

LTLRC 

79 

1  9 

0.2 

-4 . 7 

-8.8 

1  022.4 

4 . 9 

1*12 

LTLRC 

79 

110 

0  .  9 

-3 . 9 

-6 . 7 

1022.2 

7.6 

1*12 

LTLRC 

79 

111 

0  .  9 

-5.1 

-11.0 

1017.9 

14.  1 

1*12 

LTLRC 

79 

112 

0  .  9 

-2.4 

-6.8 

1  006.6 

10.9 

111 

LTLRC 

79 

113 

0.9 

-1.2 

-4 . 3 

998.8 

11.3 

112 

LTLRC 

79 

114 

0  .  0 

-7  .  1 

-15.3 

1018.2 

17.8 

1*12 

LTLRC 

79 

115 

0  .6 

-4.7 

-12.4 

1020.5 

9.0 

U2 

LTLRC 

79 

116 

1.0 

1.7 

-3.6 

1017.7 

3.2 

1*12 

LTLRC 

79 

117 

1.0 

7 . 0 

4.1 

1017.0 

2.5 

U12 

LTLRC 

79 

118 

0  .  9 

4.0 

-1.3 

1017.2 

12.0 

1*12 

LTLRC 

79 

119 

1  .  0 

6.1 

4.2 

1004.8 

5.3 

1*12 

LTLRC 

79 

120 

1.0 

7.6 

6  .  9 

993 . 1 

16.7 

1*12 

LTLRC 

79 

121 

0.7 

1.0 

-5.5 

1000.1 

19.9 

1*12 

LTLRC 

79 

122 

0.5 

0.2 

-3.3 

1006  .  1 

7.6 

U2 

LTLRC 

79 

123 

1.0 

2.2 

-0.2 

993.9 

17.8 

1*12 

LTLRC 

79 

124 

0  .  0 

-2.8 

-11.5 

1003.7 

21 . 1 

U2 

LTLRC 

79 

125 

0.7 

-1.9 

-7  .  0 

1009.1 

7  .  9 

1*12 

LTLRC 

79 

126 

1  .  0 

0  . 

-2.4 

1004.7 

10.9 

U2 

LTLRC 

79 

127 

0  .  9 

-0 . 7 

-3.8 

1001.3 

16.4 

"U2 

LTLRC 

79 

128 

0.2 

-3.1 

-9.6 

1010.3 

14.1 

1*12 

LTLRC 

79 

129 

0.6 

-3.7 

-9.7 

1013.4 

7.9 

W2 

LTLRC 

79 

130 

1.0 

0.1 

-4.6 

1010.5 

7  .  9 

W2 

LTLRC 

79 

131 

0.2 

-3.6 

-11.5 

1015.5 

19.2 

W2 

LTLRC 

79 

2  1 

0.5 

-5.3 

-11.6 

1018.1 

6.3 

1*12 

LTLRC 

79 

2  2 

1  .  0 

0.1 

-4.4 

1013.5 

3.9 

012 

LTLRC 

79 

2  3 

1.0 

3.2 

0  . 1 

1012.1 

4.6 

012 

LTLRC 

79 

2  4 

0.8 

0.3 

-5.2 

1017.0 

12.0 

1012 

LTl  RC 

79 

2  5 

1.0 

-2.2 

-10.5 

1016.6 

14.6 

012 

LTLRC 

79 

2  6 

1.0 

-0.8 

-5.3 

1004.8 

11.6 

012 

LTLRC 

79 

2  7 

1.0 

I.  1 

-2.6 

10  06.7 

8.6 

012 

LTLRC 

79 

2  8 

0.9 

0.4 

-3.5 

1012.2 

14.6 

.22 

55. 

65.5 

7.2 

3.0 

1. 

.22 

3.8 

69. 

iKttH 

7 .2 

3.0 

1  . 

.21 

0.0 

4.9 

63. 

.007 

7 . 3 

3.0 

0  .  0 

600  . 

1. 

.21 

0.0 

6  . 

57  . 

.007 

7.4 

3.0 

0  .  0 

600  . 

1  . 

0. 

C5 


142 

ITLRC 

79 

2  9 

0.2 

-8.6 

-16 . 3 

1022.8 

11.1 

142 

LTLRC 

79 

210 

0.2 

-4.2 

-9  .  9 

1018.9 

2.1 

142 

LTLRC 

79 

211 

0 . 5 

1.8 

-2.9 

1016.2 

4.2 

142 

LTLRC 

79 

212 

1.0 

8.8 

3.9 

1016.3 

7.4 

142 

LTLRC 

79 

21  J 

1.0 

1.6 

-3.3 

1015.4 

7.9 

142 

LTLRC 

79 

214 

0.7 

6.5 

4 . 0 

1  006.7 

1  .  9 

142 

LTLRC 

79 

215 

0.9 

17.6 

11.5 

1007.5 

18.5 

142 

LTLRC 

79 

216 

0.7 

1 .  0 

-6.5 

1019.8 

20.8 

142 

LTLRC 

79 

217 

1.0 

-3.0 

-8 . 5 

1020.4 

19.2 

142 

LTLRC 

79 

218 

0.7 

-0.5 

-5.4 

1020.2 

10.6 

142 

LTLRC 

79 

219 

0.2 

1 . 3 

-3.5 

1022.3 

2.3 

142 

LTLRC 

79 

220 

0.8 

3.9 

0.3 

1013.7 

4.9 

142 

LTLRC 

79 

221 

1 . 0 

10.6 

6  .  4 

1013.8 

3.9 

142 

LTLRC 

79 

222 

1  .  0 

15.7 

13.8 

1010.9 

9.7 

142 

LTLRC 

79 

22  J 

1  .  0 

17.7 

15.6 

1010.3 

4.4 

142 

LTLRC 

79 

224 

1  .  0 

7.8 

5.7 

1  0  0  7  .4 

14.8 

142 

LTLRC 

79 

225 

0.6 

4.0 
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